COMPUTING INVARIANTS OF SIMPLICIAL MANIFOLDS 



MICHAEL JOSWIG 

Abstract. This is a survey of known algorithms in algebraic topology with a focus on finite 
simplicial complexes and, in particular, simplicial manifolds. Wherever possible an elementary 
approach is chosen. This way the text may also serve as a condensed but very basic introduction 
to the algebraic topology of simplicial manifolds. 

This text will appear as a chapter in the forthcoming book "Triangulated Manifolds with Few 
Vertices" by Frank H. Lutz. 



The purpose of this chapter is to survey what is known about algorithms for the computation 
of algebraic invariants of topological spaces. Primarily, we use finite simplicial complexes as our 
model of topological spaces; for a discussion of different views see Sectional 

On the way we give explicit definitions or constructions of all invariants presented. Note that we 
did not try to phrase all the results in their greatest generality. Similarly, we focus on invariants 
for which actual implementations exist. The reader is referred to Bredon's monograph [2] for the 
wider perspective. 

For a related survey see Vegter (33] . 

1. Homology 

1.1. Definitions. Let A be a finite simplicial complex, and let A^ be the set of its /^-dimensional 
faces. Assume that the set of vertices is [n]. This way each fe-dimensional face, or k-face for short, 
AW has a unique representation {uo, i>i, • • • , v k } as an ordered set with vi < Uj+i, for which we 
simply write vqV\ ■ ■ ■ v k . 

Fixing a commutative ring R with unit element, which we write as 1, allows to define the boundary 
map 

k 

d(v vi ■ ■ ■ v k ) = ^(-1)^0 • ■ • Vi-iv i+ i ■•■v k . 

i=0 

That is to say the boundary of a k-face is a formal linear combinations of (k — l)-faces with 
coefficients in R. The set of all such linear combinations is denoted by C k _i(A; R), and its elements 
are called (k — 1)- chains. The fc-dimensional simplicial boundary operator with coefficients in R is 
the unique i?-linear extension of the boundary map to the i?-module of /c-chains C k (A;R); it is 
denoted by d k . 

The image of d k+ ±, denoted by B k (A; R), and the kernel of d k , denoted by Z k (A; R), are both R- 
modules. Their elements are called k-boundaries and k-cycles, respectively. A direct computation 
yields that d k (d k+ \(a)) = or, equivalently, B k (A;R) < Z k (A;R). The quotient 

H k (A;R) = Z k (A;R)/B k (A;R) 

is the k-th homology module of A with coefficients in R. 
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There is a certain freedom of choice in defining the map 8q. Usually, we take 8q = 0. Hence 
Hq(A; R) = R if and only if A is connected. We call a connected complex R-acyclic if all homology 
modules H k (A;R) vanish for k > 0. Contractible complexes are -R-acyclic for arbitrary R. If, 
however, we consider the empty set as a simplex of dimension —1, then the space of (— l)-chains is 
isomorphic to R as an R- module and do is surjective. By definition d-\ =0. This leads to what is 
called reduced homology: H k (A; R) = H k (A; R) for k > and H (A; R) = H (A; R) © R. 

The dimension of A, denoted as dim A, is the maximal dimension of its faces. Clearly, for 
k > dim A, we have C k (A; R) = and hence also H k (A; R) = 0. 

There are several other homology theories. Among the most important are singular and cellular 
homology. All these theories coincide for "reasonable" spaces including CW-complexes and, in 
particular, simplicial complexes; see [21 Chapter IV]. 

Not all coefficient domains are equally important. The Universal Coefficient Theorem (2J V.7.5] 
implies that knowing the integral homology of a space, that is, the case where R = Z, determines 
the homology with many other coefficient domains, including the rationals and all finite fields. Con- 
versely, the integral homology is determined by the rational homology together with the homology 
over all finite fields of prime order. Usually, we write H k (A) for H k (A;Z) and likewise for other 
Z- modules. 

1.2. Algorithms and Complexity. It is fairly easy to actually compute the homology of a space 
if the coefficient domain is a field. In this case the boundary operator is a linear map. Its kernel 
and image are vector spaces. Hence the homology modules are vector spaces whose dimensions 
can be obtained as the differences of co-rank and rank, respectively, of two subsequent boundary 
matrices. Algorithmically, the rank (or co-rank) of a matrix can be obtained by a sequence of Gauss 
elimination steps. 

We denote the number of k- faces of A by /fc(A) or simply f k . It follows that for a field F it 
requires 0(max(/fc_|_i, f k , /fc_i) 3 ) arithmetic operations in F to compute H k (A; F) provided that 
the boundary matrices d k +i and d k are given. This latter condition can be read as a warning: 
Typically a finite simplicial complex is given as the set of faces which are maximal with respect to 
inclusion. Even if their number is small (e.g., one for a simplex) the lower dimensional boundary 
matrices can still be huge: The trivial upper bound f k < m is tight. 

We define the size of A to be the number of incidences between vertices and facets. In particular, 
the size is bounded from above by (dim A + l)m, where m is the number of facets of A. Equality 
holds if and only if A is pure. This notion corresponds to an encoding of A as an array of its facets, 
where each facet is represented as a list of its vertices. 

Each integral homology module of A is a finitely generated abelian group, that is to say, a direct 
product of a finitely generated free abelian group and a product of finitely many cyclic groups of 
prime power order. The subgroup of elements of finite order is called the torsion subgroup. We 
will see below that every finitely generated abelian group, in fact, arises as a homology module of 
some 2-dimensional simplicial complex. 

On a conceptual level the computation of integral homology differs only by a little from the 
computation of rational homology. There are two reasons for this: Firstly, Q is the quotient field 
of Z. Secondly, Z is a Euclidean domain. The latter property comes in handy when one tries to 
perform a Gauss elimination over Z. Either the pivoting element is a unit, that is 1 or —1, then 
the corresponding rational basis transformations are integral, and the situation is exactly the same 
as over Q. Or, if the pivoting element is not a unit, Euclid's algorithm can be used to determine 
an integral transformation which simplifies the matrix. In this way, after finitely many steps, each 
boundary matrix d k can be transformed into diagonal shape, the Smith- Normal- Form of d k . It is 
then straightforward to determine the structure of H k (A) from the Smith-Normal- Forms of d k+ \ 
and d k . 
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The existence and also the construction of the Smith-Normal-Form of an integral matrix is a 
classical result due to H. J. S. Smith (22]; for a modern account see Munkres [HSl §11]. 

In terms of complexity there is a fundamental difference between Gauss elimination over Q and 
the sketched type of Euclidean-Gaussian elimination over Z. In the rational case the computation 
time is bounded by a polynomial (in the size of the input matrix) since the growth of the coefficients 
can be controlled. This has been observed by Edmonds jllj : see Schrijver |36l §3.3]. However, for 
the analogous operation over Z neither the size of the coefficients nor the number of arithmetic 
operations in Z is polynomially bounded. 

Relying on modular arithmetic, though, Kannan and Bachem [21] gave a first polynomial time 
Smith-Normal- Form algorithm, which was later improved by Iliopoulos |22]) Storjohann |41| . and 
others. 

While the modular approach is valid for matrices with arbitrary integer coefficients, simplicial 
boundary matrices have entries 1, —1, and only. That is to say, in an arbitrary simplicial 
boundary matrix it is always possible to perform at least a few Gauss elimination steps. Moreover, 
a typical boundary matrix is very sparse. If the matrix stays sparse during the elimination and if, 
additionally, one does not run out of unit coefficients too early (such that it is possible to continue 
with elimination steps) an elimination based Smith-Normal-Form algorithm can be superior to the 
more sophisticated methods. This is why in practical applications elimination algorithms are often 
preferred; see the survey of Dumas, Heckenbach, Saunders, and Welker ^3] and also the Section Tl. 41 
on discrete Morse theory. 

Donald and Chang [Hj proved that, asymptotically, for random simplicial complexes with many 
vertices compared to the dimension d the integral homology can be computed in 0(d/^ ax ) time 
with high probability. Here / max = max{/ , . . . , f d }. 

1.3. Examples. The homology of a 1-dimensional simplicial complex, that is, a graph, is very easy 
to understand. For a connected graph T with n vertices and e edges we immediately obtain H\(T) = 
Z e_n+1 : Consider a spanning tree, which has re — 1 edges; each one of the e — (re — 1) non-tree edges 
then corresponds to a different free generator. However, for each finitely generated abelian group G 
there is a 2-dimensional simplicial complex A such that H\(A) = G: There is a family of connected 
2-dimensional simplicial complexes C(k), for k > 2, with f(C(k)) = (1 + 3(k + 1), 3 + 12k, 9k) such 
that H 2 (C(k)) = and Hi{C{k)) ^ Z/k. From Figured which displays C(3), it should be easy to 
deduce how to define C{k) for arbitrary k; the complex C(2) is homeomorphic to the real projective 
plane RP 2 and is the only manifold in the series. The non-trivial homology cycle is carried by the 
induced subcomplex with vertices {a,b,c}. Clearly, C(k) is Q-acyclic. A wedge A V A' of two 
simplicial complexes A, A' is a complex arising by identifying one vertex in A with one vertex 
in A' in the disjoint union A U A'. It is immediate that H k (A V A'; R) = H k (A; R) H k (A'; R) for 
k > 1. By successively wedging complexes C{k) and (triangulated) cylinders S 1 x [0, 1] we obtain 
each finitely generated abelian group as the first integral homology group of a pure 2-dimensional 
simplicial complex. 

1.4. Elimination Strategies and Discrete Morse Theory. Classical Morse theory allows to 
bound the size of the homology modules of a smooth manifold by counting critical points (with 
respect to some smooth function) in a cell decomposition; see Milnor |29l §1.5]. Forman |13l I14j 
developed a Morse theory for abstract simplicial complexes (and, more generally, for arbitrary finite 
CW- complexes), which we briefly sketch now. 

Let A be a finite simplicial complex, and let p : A — > N be an arbitrary function. For each 
simplex a 6 A we define the critical defect 

<W = # {P S A | p < a, n(p) > n(a)} + # {r € A | r > a, fi(r) < p{a)} , 
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Figure 1. Complex C(3): All unlabeled vertices are distinct. The three special 
vertices a, b, c occur three times each in the boundary of the 2-cell suggested by the 
picture. 



where by p < a we denote that p is a maximal face of a. Now a function p : A — » N is a discrete 
Morse function of A if for each simplex the critical defect is at most one. A simplex of critical 
defect zero is called critical; otherwise it is regular. Discrete Morse functions always exist: The 
function which maps each simplex to its dimension is a Morse function for which each simplex is 
critical. 

Fix a discrete Morse function p. Then, by definition, for each regular simplex a there is either a 
unique regular face p < a with p(p) > p(a) or a unique regular face r > a with p(r) < p(a). That 
is to say, each discrete Morse function gives rise to a matching in the Hasse diagram of the face 
lattice of A, the Morse matching of p. Each Morse matching [p] satisfies the following property: 
Direct all edges of the Hasse diagram of A consistently, say, with increasing dimension, and reverse 
the orientation for each edge in [p\; then there is no oriented cycle. Conversely, each matching 
satisfying this acyclicity condition is associated with a discrete Morse function. For an example of 
a Morse matching in the triangulation RPg of the real projective plane see Figures El El an d IU 

Each edge o~<t in the Hasse diagram of the face lattice of A, for r £ A^ (and thus a G A^ fc_1 )), 
corresponds to a ±1 entry in the boundary matrix dk and hence can be used to eliminate all other 
entries in the row of r and in the column of a. Now, a subset E of the edges of the Hasse diagram 
corresponds to a sequence of elimination steps which can be performed (in an arbitrary order) 
without interfering with each other if and only if E is a Morse matching. That is to say, a globally 
optimal elimination strategy corresponds to a discrete Morse function with as few critical faces as 
possible. 

Figure |1] shows the Morse matching [p] of MPg from Figures [21 and |3] as a set of entries in the 
boundary matrices. Interpreting each elimination step geometrically leads to a transformation 
from a triangulation to a more general finite CW-complex. In this particular case, we obtain a 
cell decomposition of RP 2 into three cells or dimensions 0, 1,2, respectively. The unique one cell 
corresponds to a (real) projective line. This cell decomposition is minimal since a CW-complex 
with only two cells (one of which must be a point, by definition) necessarily is a sphere. Further 
this argument implies that the Morse matching p of RPg (which leads to this minimal cell decom- 
position) is an optimal Morse matching. Note that Forman more generally proved |14l Corollary 
1.7] that a combinatorial manifold without boundary which admits a Morse matching with only 
two critical cells is a shellable sphere. In this sense optimal Morse matchings generalize shelling 
orders to arbitrary simplicial complexes. 
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Figure 2. Triangulation MPg of the projective plane (obtained from identifying an- 
tipodal points in the boundary of the regular icosahedron) with a Morse matching \p] 
(indicated by the arrows). This triangulation is minimal with /-vector (6,15,10). 
The regular faces are marked with a circle or a square, depending whether they are 
matched to higher or lower dimensional face. There are three critical faces: 4, 45, 
and 012. It is clear that fi cannot be extended: The critical triangle 012 cannot 
be matched since all its edges are already matched; the critical vertex 4 cannot be 
matched to the critical edge 45 since this would violate the acyclicity condition. The 
subsequent Figures 01 and 0] display the same information in different ways. 




FIGURE 3. The Morse matching [fj] in the Hasse diagram of the face lattice of KPg. 
The bold edges correspond to the arrows in Figure [51 The critical faces are un- 
matched. 



While optimal discrete Morse functions are important, it is not easy to find good discrete Morse 
functions algorithmically. 

Theorem 1. (Egecioglu and Gonzalez For the class of finite pure 2- dimensional simplicial 
complexes which are embeddable into M 3 it is NP-complete to decide whether there is a Morse 
matching with at most c critical 2-faces. 
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Figure 4. Boundary matrices d 2 and d\ of MPg. The matrix d\ is transposed in 
order to visualize how the two matrices fit together. The boxed entries in the two 
matrices correspond to the arrows in Figure [21 and the matching edges in Figure |3 
Using the boxed entries as pivots in an arbitrary order yields an elimination strategy 
for the computation of the homology modules. 

In fact, Egecioglu and Gonzalez even proved a rather strong non-approximability result, which 
we omit. Notice that, for the result above and the corollary below, the number c is part of the 
input. 

One can show that an optimal discrete Morse function of any connected simplicial complex 
necessarily has a unique critical vertex. This can be used to obtain the following result |23j : see 
also Lewiner |27j . 

Corollary 2. For the class of finite pure 2 -dimensional simplicial complexes which are embeddable 
into M 3 it is NP '-complete to decide whether there is a Morse matching with at most c critical faces 
of arbitrary dimension. 

For a more stringent interpretation of discrete Morse theory from the viewpoint of shellings see 
Chari [3]- The set of all optimal Morse functions of a given complex may have an interesting 
structure, although in most cases it is very difficult to obtain; see |S] for a discussion of special 
cases. 

In practice a greedy choice of an elimination strategy /discrete Morse function often leads to a 
considerable reduction of the input to a homology computation which is then passed on to some 
Smith-Normal-Form algorithm. This way it is possible to compute the homology of rather large 
simplicial complexes; e.g., see [To] . 

1.5. Bases Transformations. In order to locate subcomplexes which carry the homology cycles 
it is necessary to keep track of all the transformations of bases during an elimination process. Even 
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if the boundary matrix stays sparse, the accumulated bases transformation matrices almost never 
do. If we want to avoid to run out of memory during the computation it is therefore useful to 
store the transformations as sequences of elementary transformations, which in turn are stored in 
compressed form. Such techniques are also common in numerical linear algebra. 

2. COHOMOLOGY 

Transposing each boundary matrix leads to coboundary maps, cochains, cocycles, and, finally, 
cohomology modules, which contain the same information as the homology modules, but in a dual 
form. Everything concerning the computation of homology modules can directly be translated to 
cohomology. Here we are interested in the cup product, which has a less natural description in 
terms of homology. 

2.1. Definitions. Let again A be a finite simplicial complex, and let R be a commutative ring. 
The -R-module of /c-dimensional cochains is the set C k (A;R) of /2-linear maps from C k (A) to R. 
The evaluation of / € C k (A;R) at a € C k (A;R) is written as (f,a). For a = v ■ ■ ■ v k € A (fc) 
let a* be the /c-cochain with maps a to 1 and all other fc-simplices to 0, and for each face r = 
vq ■ ■ ■ ViWi + i ■ ■ ■ v k € A( fc+1 ) let e(<r, t) = (— \) t+l . The latter value is the coefficient of a in the 
boundary of r. We define the coboundary of a* as 

t7><j 

Via the canonical identification a <-> a* we have that the i?-linear extension 5 k : C k (A;R) —> 
C k+l {A; R) of 5 is the transpose of the map 

Again, the image of 5 k ~ 1 , denoted by B k (A;R), and the kernel of 5 k , denoted by Z k (A;R), 
are both iZ-modules. Their elements are called k-coboundaries and k-cocycles, respectively. The 
quotient 

H k (A;R) = Z k (A;R)/B k (A;R) 

is the A;-th cohomology module of A with coefficients in R. 

One reason to look at cohomology rather than homology is that cohomology comes equipped 
with a natural ii-bilinear function, the cup product U : W{A- R) x H k {A; R) -> H i+k {A; R) which 
is defined on the cochain level by 

(/ U g, v ■ ■ ■ v i+k ) = (-l) lk (f,v ---Vi) (g,Vi---v i+ k), 

where / £ C^A-R) and g G C k {A;R). This turns H*(A;R) = Q) k H k (A;R) into an associative 
-R-algebra, the cohomology algebra of A. Note that the cup product is not commutative, in general. 
Instead we have 

(1) 9 Uf = (-iy k fUg. 

There is another interesting product, the cap product n : H l (A;R) x H k (A;R) — > H k _i(A; R), 
which can be defined, again on the level of cochains and chains, by 

/ n v ■ ■ ■ v k = (-!)*(*-<) (/, VQ ... Vi ) Vi ... Ufcj 

for i < k. For details see [21 §VI.5]. 
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2.2. Algorithms. For the integral homology modules it was obvious what the output of a com- 
putation should be: Each finitely generated abelian group has a canonical finite representation by 
its free rank and its torsion coefficients. The goal for computing the cup product is somewhat less 
obvious: Computing the cohomology modules is the same as for homology. If we keep track of all 
basis transformations we also obtain a set of cochains which generates H*(A;R) as an R- module. 
This in turn can be used to compute a multiplication table for the cup product; see Peitgen |34| . 

However, there is no canonical choice for a generating system. This has the serious disadvantage 
that there does not seem to be a direct way to decide whether two cohomology rings are isomorphic 
or not. 

2.3. Duality of Manifolds. Let M be a closed d-dimensional simplicial manifold. Since each 
codimension-l-face is contained in exactly two facets we have H d (M; Z/2) = Z/2. We denote the 
generator of H d (M; Z/2) by [M] 2 . The manifold M is oriented if H d (M) ^ Z. If M is not oriented 
then H d (M) = 0. 

The following theorem comes in many guises: 

Theorem 3. (Poincare Duality; see PI VI. 8. 4]) The map 

H k (M;Z/2) -» H d _ k (M;Z/2) : / » fn [M} 2 

is an isomorphism. If, additionally, M is oriented, and if [M] is a generator of H d (M), then the 
map 

H k (M) — ► H d _ k (M) :/•->/ D [M] 

is an isomorphism. 

2.4. Stiefel- Whitney Classes. Let M be a closed d-dimensional smooth combinatorial manifold. 
The fc-th Stiefel- Whitney class of M is a certain (characteristic) cohomology class u k G H k (M; Z/2). 
For a precise definition see ;2 4 §VI.17]. The k-th. Stiefel- Whitney class forms the essential obstruction 
to the existence of d—k+l linearly independent cross sections in the tangent bundle of M\ see Milnor 
and Stasheff |31[ Proposition 4.4]. In the sequel we give an elementary combinatorial description 
of io k due to Goldstein and Turner JH]; see also Halperin and Toledo |2l"] . 

For a pair of faces a = vq ■ ■ ■ v k and r with a < r define the i-th gap as 

Gi(a, t) = {w G r | Vi < w < v i+1 } , 

where < i < k. Additionally, let G-i(a, r) = {w G r | w < vq} and G k = {w G r | v k < w}. 
Now (a, t) is a regular pair if for all i odd Gi(a,r) = 0. 

Let u>k be the ^-dimensional chain with coefficients in Z/2 which is the sum of all k- faces a, for 
which the number of /-faces r, such that I > k and (a, r) is a regular pair, is odd. It is easy to 
verify that uj k is a homology cycle mod 2. We identify uj k with its Z/2-homology class. 

Since for each A;-face a the pair (cr, a) is regular, it is obvious that io d = [Af]a 7^ 0. 

Theorem 4. (Goldstein and Turner |19j ) The Z/2-homology class ofcuk is the image of the (d—k)- 
th Stiefel- Whitney class u> d ~ k under the Poincare Duality isomorphism. 

2.5. The Intersection Form of a 4-Manifold. Let M be a closed oriented combinatorial 4- 
manifold. This implies that M carries a unique smooth structure. 

Since M is oriented we have H^(M) = Z. We fix a generator [M] G H^(M), which is called an 
orientation class. In view of equation Q the map 

Q M : H 2 (M) x H 2 (M) -» Z : (/, g) ^ (f U g, [M] ) 

is a symmetric bilinear form, the intersection form of M. If / or g is of finite order, then Qm(/ 5 <?) 
( ). Therefore, by choosing a basis, we can identify Qm with an integral (r x r)-matrix M, the Gram 
matrix of Qm, where r is the free rank of H 2 (M). It can be shown that Qm is unimodular, that is, 
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det Qm = ±1- The form Qm is even if all diagonal elements of M are even and odd otherwise. Note 
that this property does not depend on the basis chosen. Over the reals Qm can be transformed 
into a diagonal matrix D with entries ±1. The number sign(Qj\/) = (number of 1-entries) — 
(number of —1-entries) in D is the signature of Qm- If Qm is indefinite, that is, |sign(QM)| 7^ r, 
then up to an integral basis transformation Qm is uniquely determined by the rank, the parity, and 
the signature; see Milnor and Husemoller |301 Theorem II. 5. 3]. The situation is more complicated 
if Qm is a (positive or negative) definite form; see [3UJ §11.6]. 

The Wu formula [201 Proposition 1.4.18] implies that the intersection form of M is even if and 
only if the second Stiefel- Whitney class vanishes. 

Whitehead proved that the homotopy type of a simply connected closed 4-manifold is determined 
by its intersection form; see [301 Theorem V.1.5]. However, this has been improved to a considerably 
stronger famous result. 

Theorem 5. (Freedman ^]) Two simply connected closed oriented combinatorial ^-manifolds are 
homeomorphic if and only if their intersection forms are equivalent. 

If we additionally assume that M is simply connected, then H\{M) = and hence, by the 
Universal Coefficient Theorem V.7.5], H (M) = 0. The Poincare Duality Theorem |31 now gives 
that Hs(M) = 0. Moreover, H2(M) = H 2 (M) is torsion-free and thus the intersection form Qm 
completely determines the cohomology ring of M. 

The intersection form can also be defined in terms of the cap product (see (21 Corollary VI. 5. 3]) 

(2) (fUg,[M}) = (f,gn[M]), 

where f,g£ H 2 (M). Recall that the map g g H [M] is the Poincare Duality isomorphism. 
For more information on the topology of 4-manifolds see Kirby |25j . Gompf and Stipsicz |2()j . 

2.6. An Example. The following is the list of facets of the unique minimal triangulation CPg of 
the complex projective plane with /-vector (9,36,84,90,36); see Kiihnel and Banchoff |26| : 

01234 01236 01247 01267 01345 01356 01457 01568 01578 

01678 02348 02367 02378 02457 02458 02578 03456 03468 

03678 04568 12348 12356 12358 12467 12468 12568 13457 

13478 13578 14678 23567 23578 24567 24568 34567 34678 

We sketch how to prove that this is, indeed, a triangulation of CP 2 . As a first step we have 
to prove that this is actually a combinatorial 4-manifold. To this end one has to verify that the 
link of each /c-face is a (3 — /c)-sphere, the crucial step being k = 0. Recognizing the 3-sphere is 
algorithmically possible but the methods proposed by Rubinstein, Thompson [12]) and others do 
not seem to be feasible in practice. Fortunately, however, for small cases, such as this example, the 
flip-heuristics by Bjorner and Lutz solves the problem easily. 

Since /^(CPg) = (3) the triangulation is 3-neighborly and hence the space is simply connected. 
In view of Theorem it then suffices to determine the intersection form. 

Using elimination techniques we obtain that 

H°(CPl) ^ H 2 (CPl) ^ F^CP 2 ,) ^ Z 

and all other cohomology groups vanish. Additionally, by keeping track of the bases transforma- 
tions, one can obtain a generator for H 2 (CPg), for instance, b = 125* + 128* + 137* + 138* + 146* + 
148* - 235* + 246* + 256* - 268* - 347* - 357* - 358* + 467* - 478* + 567*. Instead of computing 
the cup product b U b directly, we first compute a generator for #2 (CPg) as 

b n [CP^ = 567 - 568 + 578 - 678. 
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From this we obtain Q CP 2(b, b) = (b, b n [CP|]) = 1. In particular, 6 U b generates i? 4 (CP|). We 
conclude that the intersection form Q^p2 has a 1 x 1-Gram matrix (1) and that 

H*(CPl) ^Z[b]/(b 3 ). 

The homology 2-cycle bd [CPg] corresponds to the boundary of the non-face 3-simplex 5678, that 
is, #2 (CPg) is generated by an embedded 2-sphere. Note that this is typical in the sense that every 
element of the second homology group of a closed, oriented, smooth 4-manifold is generated by an 
embedded surface; see [201 Proposition 1.2.3 and Chapter 2]. In our case, 6 Pi [CPg] corresponds to 
a complex projective line. 

Now we turn to the computation of the Stiefel- Whitney homology classes of CPg. By |2f)| 
Proposition 1.4.18] the result Q CF 2(b,b) ^ above already implies that w 2 does not vanish. A 
direct computation can be sketched as follows. Enumerating all regular pairs yields that uj\ and LV3 
vanish, uj = + 5 + 8, and that lo 2 = 012 + 018 + 023 + 027 + 028 + 034 + 045 + 056 + 057 + 058 + 
067 + 078 + 126 + 127 + 128 + 146 + 147 + 238 + 246 + 247 + 257 + 258 + 348 + 357 + 358 + 378 + 
456 + 457 + 458 + 467 + 567 + 568 + 578 + 678. Reducing modulo boundaries simplifies the result 
to u>2 = 567 + 568 + 578 + 678 (and uiq is the Z/2-homology class of a point). 

3. Fundamental Group 

For an arbitrary topological space T and a point * € T let iti(T, *) be the set of homotopy 
classes of closed paths based at *. The concatenation of paths imposes a group multiplication on 
7Ti (T, *). This gives the fundamental group of T with base point *. For pathwise connected T the 
structure of the fundamental group does not depend on the choice of the base point. In this case 
we omit the base point in the notation. The space T is simply connected if it is pathwise connected 
and vti(T) is the trivial group. 

Now let A be a connected finite simplicial complex, and let S C A^ 1 ) be the set of edges of a fixed 
spanning tree. We define a finitely presented group E(A,S) by listing generators and relations: 
Each edge vqVi € A^ 1 ) \ S is a generator, and for each 2-face woWiw 2 G A^ 2 ) we have a relation 

(w w 1 )(w 1 w 2 )(w w 2 y 1 = 1, 

with the understanding that for WiWk € S the corresponding factor is the identity. Then E(A, S) = 
7r(A). For a combinatorial approach to fundamental groups see Seifert and Threlfall |37| Chapter 7]. 

It is easy to show that each finitely presented group arises as the fundamental group of a finite 
simplicial complex. In fact, it suffices to suitably modify the construction in Section ll.3l of a sim- 
plicial complex with a given finitely generated abelian group as its first homology; see Stillwell |40| 
§3.4.4]. Each finitely presented group even arises as the fundamental group of a 4-manifold. 

Therefore, in spite of the fact that it is simple to obtain a finite description for the fundamental 
group of A, this is often not very useful: In general, it is not possible to decide whether a given 
finitely presented group is finite or trivial. Moreover, it is, in general, impossible to decide whether 
two finitely presented groups are isomorphic. As a consequence it is impossible to decide whether 
two finite simplicial complexes are homotopy equivalent or homeomorphic. For these insolvability 
results see [113 Chapter 9]. 

However, there is an algorithm, due to Todd and Coxeter [33], which, for a given finitely presented 
group G and a subgroup U of finite index, enumerates all cosets U g with g € G. In particular, if 
G happens to be finite, then the Todd-Coxeter algorithm can be used to determine the order of G. 
This and related algorithms can be used to obtain simplified presentations in particular cases. 

The existence of coverings is tightly related to the structure of the fundamental group of the 
base space; see Chapter III] and, for a constructive approach, Rees and Soicher [35] ■ 
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The first homology group ifi(A) coincides with the quotient of 7i"i(A) by its commutator group. 
In particular, if 7Ti(A) is abelian then 7Ti(A) = Hi(A). In this case the Todd-Coxeter algorithm 
can, of course, be replaced by any Smith-Normal-Form algorithm. 

Brown [I] gave an algorithm to compute the higher homotopy groups of any simply connected 
finite simplicial complex; see also Sergeraert |381 Theorem 10.2]. 

4. Data Types and Other Combinatorial Models for Spaces 

So far we chose a finite simplicial complex as our model of a topological space. This allowed to 
access several algebraic invariants algorithmically. Two questions arise: 

(1) Is there a more efficient description? 

(2) Is there a computationally feasible model which comprises more spaces? 

The first question is easily motivated by looking at the boundary dAd of a single d-simplex: 
Since /^(SA^) = the size of the boundary matrices is not bounded by a polynomial in d, 

and hence the time required to compute the simplicial homology (via the methods sketched in this 
survey) is exponential in d. On the other hand, dA^ is homeomorphic to which has a cell 

decomposition with two cells only. Applying the same Smith-Normal-Form algorithms as for the 
simplicial theory to the cellular boundary matrices allows to compute all cellular homology modules 
of dAd in linear time. For an introduction to cellular homology see [2J §IV.10]. 

The problem here is to determine the boundary matrices. Since the characteristic maps, which 
describe how (the boundary of) a particular fc-cell is glued to the (k — l)-skeleton, can be wild, 
there is no hope for a purely combinatorial description of arbitrary cell complexes. So one solution 
could be to represent a cell complex by its cellular boundary matrices. However, this means to 
loose (quite a lot of) information: In general, it is not even possible to define the cup product from 
the cellular boundary matrices. 

There are special types of cell complexes with restricted characteristic maps, which have a useful 
combinatorial description: A regular cell decomposition of a surface can be represented, e.g., as a 
quad-edge data structure. For an overview see O'Rourke [331 §4.4.3]. For a certain generalization 
to higher dimensions see Brisson (3]. 

One other useful special class of regular cell complexes are the pseudo-simplicial complexes: each 
cell is a simplex, but any two simplices are allowed to share several faces in their boundary. For 
instance, a sphere of arbitrary dimension d can be represented by two d-simplices which share their 
whole boundary. 

Since every cell complex can be triangulated, replacing finite simplicial complexes by general 
finite cell complexes does not enlarge the class of spaces accessible. Employing techniques from 
functional programming Sergeraert suggested [3Ej to encode a space T as a function which, for 
each k, computes the fc-th homology group. His class of spaces with effective homology, e.g., is 
closed under forming loop spaces. This way it is feasible to compute the higher homotopy groups 
of simply connected finite simplicial complexes |38| Theorem 10.2]; see also 

5. Software 

A program to compute simplicial homology is a (proposed) GAP package due to Dumas et 
al. As a special feature this program offers a variety of algorithms, including modular ones. 

De Silva's PLEX [7j is a collection of MATLAB functions to be used as a research tool for 
building and studying simplicial complexes. This includes Z/2-homology computation as well as 
rational homology (with floating point arithmetic). 

SnapPea is a program for creating and studying hyperbolic 3-manifolds |45| . It focuses on metric 
properties. 
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The flip-heuristics of Bjorner and Lutz to simplify a given triangulation of a space have been 
implemented by Lutz |28j . 

polymake has originally been designed as a software package for studying convex polytopes |17l 
118) . Yet, starting with version 2.0, polymake's application TOPAZ deals with finite simplicial com- 
plexes. Besides many standard and a few non-standard constructions, it primarily supports to 
compute algebraic invariants: simplicial homology and cohomology modules (with and without 
bases, via elimination), cap-products, cup-products, Stiefel- Whitney classes, and intersection forms 
of 4-manifolds. In fact, the example computation for the complex projective plane, see Section \2. 61 
had been performed with polymake/TOPAZ. 
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